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Abstract 

We study the solid-on-solid interface model above a horizontal wall 
in three dimensional space, with an attractive interaction when the 
interface is in contact with the wall, at low temperatures. The system 
presents a sequence of layering transitions, whose levels increase with 
the temperature, before the complete wetting above a certain value of 
this quantity. 
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Consider the square lattice Z 2 . To each site x = (aq,^) of the lattice, 
an integer variable <p x > is associated which represents the height of the 
interface at this site. The system is first considered in a finite box A C Z 2 
with fixed values of the heights outside. Each interface configuration on A: 
{4> x }, x G A, denoted 0a, has an energy denned by the Hamiltonian 

# A (0 A |0) = 2J \<l> x -(f>x>\-2(J-K)J2$((f>x) + 2J\M, (1) 

where J and K are positive constants, the function 8 equals 1, when (p x = 0, 
and 0, otherwise, and |A| is the number of sites in A. The first sum is taken 
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over all nearest neighbors pairs (x,x ! ) C Z 2 , such that at least one of the 
sites belongs to A, and one takes <p x = <p x when x (jL A, the configuration <fi 
being the boundary condition, assumed to be uniformly bounded. 

In the space R 3 , the region obtained as the union of all unit cubes centered 
at the sites of the lattice A x Z, that satisfy x 3 < 0(a?i, x 2 ), is supposed to be 
occupied by fluid +, while the complementary region above it, is occupied by 
fluid — . The common boundary between these regions is a surface in R 3 , the 
microscopic interface. The region £3 < —1/2 is considered as the substrate, 
also called the wall W. 

The considered system differs from the usual SOS model by the restriction 
to non-negative height variables and the introduction of the second sum in 
the Hamiltonian, the term describing the interaction with the substrate. 

The probability of the configuration <f) A , at the inverse temperature = 
1/kT, is given by the finite volume Gibbs measure 

ti A (<f> A I 0) = Z(A, exp ( - /3H A (<f) A I $)), (2) 

where Z(A, <fr) is the partition function 

Z(A,$) = £exp ( _ f3H A (4> A I $)). (3) 

Local properties at equilibrium can be described by correlation functions 
between the heights on finite sets of sites, obtained as expectations with 
respect to the Gibbs measure. 

We next briefly describe some general results, which are an adaptation to 
our case of analogous results established by Frohlich and Pfister (ref. [Tj) for 
the semi-infinite Ising model. 

Let A C Z 2 be a rectangular box of sides parallel to the axes. Consider 
the boundary condition (p x = 0, for all x (jL A, and write Z(A, 0) for the 
corresponding partition function. The associated free energy per site, 

r w - = - lim(l//3|A|)lnZ(A,0), (4) 

A— >oo 

represents the surface tension between the medium — and the substrate W . 
This limit @ exists and < t w ~ < 2 J. One can introduce the densities 

p z = lim £ (5(<p x - z'))f\ po = lim (5(<p x ))f\ (5) 

A— >oo , — ' A— >-oo 
z'=0 
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Their connection with the surface free energy is given by the formula 



The surface tension t w+ between the fluid + and the substrate is t w+ = 0. 
In order to define the surface tension t + ~ associated to a horizontal inter- 
face between the fluids + and — we consider the ordinary SOS model, with 
boundary condition <j) x = 0. The corresponding free energy gives r + ~. With 
the above definitions, we have 



and the right hand side in (jTJ) is a monotone increasing and concave (and 
hence continuous) function of the parameter K. This follows from relation 
fl6]) where the integrand is a positive decreasing function of K. Moreover, 
when K > J equality is satisfied in (JTj). 

In the thermodynamic description of wetting, the partial wetting situation 
is characterized by the strict inequality in equation (J7|), which can occur 
only if K < J, as assumed henceforth. We must have then p > 0. The 
complete wetting situation is characterized by the equality in ([7]). If this 
occurs for some K, say K' < J, then equation (jH]) tells us that this condition 
is equivalent to p = 0. Then, both conditions, the equality and p = 0, hold 
for any value of K in the interval (K', J). 

On the other hand, we know that po = implies also that p z = 0, for any 
positive integer z. This indicates that, in the limit A — > oo, we are in the 
+ phase of the system, although we have used the zero boundary condition, 
so that the medium — cannot reach anymore the wall. This means also that 
the Gibbs state of the SOS model does not exist in this case. 

That such a situation of complete wetting is present for some values of 
the parameters does not follow, however, from the above results. Actually 
this fact, as far as we know, remains an open problem for the semi-infinite 
Ising model in 3 dimensions. For the model we are considering an answer to 
this problem has been given by Chalker [2]. 

Chalker's result. We use the following notation: 





T W+ ((3)+T + -(P)>T W -(P,K). 



(7) 



u = 2p(J-K), t = e" 4/3J . 



(8) 



If u < - ln(l - t 2 ), then p = 0. 
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Thus, for any given values of J and K , there is a temperature below which 
the interface is almost surely bound and another higher one, above which it 
is almost surely unbound and complete wetting occurs. At low temperature 
(i.e., if u > In 16), we have p > 0. 

The object of our study is to investigate the region not covered by these 
results when the temperature is low enough. As mentioned in the abstract, 
we shall prove that a sequence of layering transitions occurs before the system 
attains complete wetting. More precisely the main results can be summarized 
as follows. 

Theorem 1. Let the integer n > be given. For any e > there exists 
a value t (n, e) > such that, if the parameters t, u, satisfy < t < t (n, e) 
and 

- ln(l - t 2 ) + (2 + e)r+ 3 < u < - ln(l - t 2 ) + (2 - e)t n+2 , (9) 

then the following statements hold: (1) The free energy t w ~ is an analytic 
function of the parameters t,u. (2) There is a unique Gibbs state p n , a 
pure phase associated to the level n. (3) The density is po > 0. The second 
inequality in is not needed in the case n = 0. 

An illustration for this theorem, in the plane (K, is given in Figure 
1. From it we can see, as mentioned in the abstract, that if the parameter 
K is kept fixed, that seems natural since it depends on the properties of the 
substrate, then the value n of the level increases when the temperature is 
increased. 

Concerning this theorem, the following remarks can be made: 

(1) The analyticity of the free energy comes from the existence of a con- 
vergent cluster expansion for this system. This implies the analyticity, in a 
direct way, of some correlation functions and, in particular, of the density p . 

(2) The unicity of the Gibbs state means that the correlation functions 
converge, when A — > oo, to a limit that does not depend on the chosen 
(uniformly bounded) boundary condition (p x . Being unique and translation 
invariant this state represents a pure phase. It is associated to a level n in 
the sense that, for the typical configurations of the state, large portions of 
the interface are near to the level n. 

(3) The condition po > means that the interface remains at a finite 
distance from the wall and hence, we have partial wetting. We can see that 
the region where this condition holds is, according to the Theorem, much 
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larger, at low temperatures, than the region initially proved by Chalker. It 
comes very close to the line above which it is known that complete wetting 
occurs. 

(4) We have, t (n, e) — >■ when n — > oo or e — > 0. 

The reason why t (n, e) depends on e, satisfying remark 4, has an expla- 
nation. One may believe that the regions of uniqueness of the state extend 
in such a way that two neighboring regions, say those corresponding to the 
levels n and n + 1, will have a common boundary where the two states ji n 
and fi n+ i coexist. 




FIGURE 1. The analyticity regions of Theorem 1. 

At this boundary there will be a first order phase transition, since the 
two Gibbs states are different. The curve of coexistence does not exactly 
coincide with the curve u — — ln(l — t 2 ) + 2t ra+3 . Theorem 1 says that it is 
however very near to it, if the temperature is sufficiently low. 

Let us formulate in the following statement the kind of theorem that 
we expect. We think that such a statement could be proved using, as for 
Theorem 1, an extension of the Pirogov-Sinai theory. 

Statement. For each given integer n > 0, there exists t (n) > and a 
continuous function u = ip n+ i{t) on the interval < t < t (n), such that 
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the statements of Theorem 1 hold, for t in this interval, in the region where 
ip n+ i{t) < u < i/) n+2 (t) and for n = 0, in the region ipi(t) < u. When 
u = ip n+ i(t) the two Gibbs states, \i n and (i n +i, coexist. 

The existence of a sequence of layering transitions has been proved for a 
related model, known as the SOS model with an external magnetic field. See 
the works by Dinaburg, Mazel [3], Cesi, Martinelli jl] and Lebowitz, Mazel 
[5]. This model has the same set of configurations as the model considered 
here, but a different energy: The second term in (PQ) has to be replaced by the 
term +h Y^xeA (px to obtain the Hamiltonian of the model with an external 
magnetic field. The method followed for the proof of the Theorem is essen- 
tially analogous to the method developed for the study of that model. The 
most important difference between the two systems concerns the restricted 
ensembles and the computation of the associated free energies. 

Concerning the proof of Theorem 1 (paper in preparation) let us say that, 
for an interesting class of systems, among which our model is included, one 
needs some extension of the Pirogov-Sinai theory of phase transitions (see 
ref. [6]). In such an extension certain states, called the restricted ensembles, 
play the role of the ground states in the usual theory. They can be defined 
as a Gibbs probability measure on certain subsets of configurations. In the 
present case one considers, for each n = 0, 1, 2, . . ., subsets of configurations 
which are in some sense near to the constant configurations <p x = n. 

Namely, we consider the set C£ es (A, n) of the microscopic interfaces, with 
boundary at height <f) x = n, and whose Dobrushin walls have, all of them, 
horizontal projections with diameter less than 3k + 3 (these walls are the 
maximally connected sets of vertical plaquettes of the interface). The Gibbs 
measure defined on the subset C£ es (A,n) is the restricted ensemble corre- 
sponding to the level n. The associated free energies per unit area 

/ fc (n) = - lim (1//?|A|) In £ exp(-/3F(0 A |n)) (10) 

A eq es (A,n) 

can be computed, with the help of cluster expansions (see, for instance, 
ref. [7]), as a convergent power series in the variable t. Then one is able to 
study the phase diagram of the restricted ensembles. The restricted ensemble 
at level n is said to be dominant, or stable, for some given values of the 
parameters u and t, if fk{n) = min n / fk(n'). We then have: 

Proposition 1. Let the integer n > be given and choose k > 1. Let 
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a, b > be two real numbers. Let < t < t\(k) — (3k + 3) 4 . If 

- ln(l - t 2 ) + (2 + a)t n+3 < u < - ln(l - t 2 ) + (2 - 6)t n+2 , (11) 
t/ien 7 we have 

fk(n) < f k (h) - at 3n+3 + 0(t 3n+4 ), foranyh>n + l, (12) 
/fcW < /fc(/i) -&t 3n + 0(t 3n+1 ), for anyO<h<n-l. (13) 

We notice that A; > n is the useful case in the proof of Theorem 1, and 
that the remainders in inequalities ( Tl2l and ( TT3]) can be bounded uniformly in 
h. Then the proof of Theorem 1 consists in showing that the phase diagram 
of the pure phases at low temperature is close to the phase diagram of the 
dominant restricted ensembles. 
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